UNCLASSIFIED 


DEFENSE  DOCUMENTATION  CENTER 

FOR 

SCIENTIFIC  AND  TECHNICAL  INFORMATION 

CAMERON  STATION.  ALEXANDRIA,  VIRGINIA 


UNCLASSIFIED 


NOTICE:  When  government  or  other  drawings,  speci¬ 
fications  or  other  data  are  used  for  any  purpose 
other  than  in  connection  with  a  definitely  related 
government  procurement  operation,  the  U.  S. 
Government  thereby  incurs  no  responsibility,  nor  any 
obligation  whatsoever;  and  the  fact  that  the  Govern¬ 
ment  may  have  formulated,  furnished,  or  in  any  way 
supplied  the  Baid  drawings,  specifications,  or  other 
data  is  not  to  be  regarded  by  implication  or  other¬ 
wise  as  in  any  manner  licensing  the  holder  or  any 
other  person  or  corporation,  or  conveying  any  rights 
or  permission  to  manufacture,  use  or  sell  any 
patented  invention  that  may  in  any  way  be  related 
thereto. 


THIS  DOCUMENT  IS  BEST 
QUALITY  AVAILABLE.  THE  COPY 
FURNISHED  TO  DTIC  CONTAINED 
A  SIGNIFICANT  NUMBER  OF 
PAGES  WHICH  DO  NOT 
REPRODUCE  LEGIBLY. 


434192 


Columbia  ^nibertfitp 

intbfCilpotHetoJorh 

DEPARTMENT  OF  CIVIL  ENGINEERING 
AND  ENGINEERING  MECHANICS 


FORCED  VIBRATIONS  OF  A  VISCO-ELASTIC  CYLINDER 
CASE-BONDED  TO  A  THIN  ELASTIC  SHELL 


by 

L.  A.  Henry  and  A.  M.  Freudenthal 


Office  of  Naval  Research 
Project  NR  064-446 
Contract  Nonr  266(78) 
Technical  Report  No.  22 


January  1964 


Reproduction  in  whole  or  in  part  is  permitted  for  any  purpose 
of  the  United  States  Government. 


Columbia  ©niberoitp 

intljtCiipotllrtogorft 

DEPARTMENT  OF  CIVIL  ENGINEERING 
AND  ENGINEERING  MECHANICS 


FORCED  VIBRATIONS  OF  A  VISCO-ELASTIC  CYLINDER 

/ 

CASE-BONDED  TO  A  THIN  ELASTIC  SHELL 

/ 


k 


by 


L.  A.  Henry  and  A.  M. 


Freudenthal, 


Office  of  Naval  Research 
Project  NR  0644446 
Contract  Nonr  266(78) 
Technical  Report  No.  22 


January  1964 


Reproduction  in  whole  or  in  part  is  permitted  for  any  purpose 
of  the  United  States  Government. 


► 


ABSTRACT 


4 

This  paper  is  concerned  with  the  forced  vibra¬ 
tions  of  case-bonded  soliu  propellant  grains. 
Special  attention  is  paid  to  the  interaction  be¬ 
tween  the  thin  elastic  shell  (case  and  the  thick- 
walled  propellant  cylinder.  The  effects  of  ma¬ 
terial  dampiny  in  the  propellant  are  explored. 

Frequency- response  functions  are  computed, 
which  can  be  used  for  analysis  in  situations  where 
the  loading  is  arbitrary  deterministic  or  random 
stationary  . 
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I .  INTRODUCTION 


Increasing  interest  in  the  dynamic  response  of  solid 
propellant  rocket  motors  has  stimulated  considerable  research 
in  the  area  of  vibrating  cylinders.  Many  papers  have  been 
written  dealing  with  the  dynamic  response  of  thick-walled 
elastic  cylinders  and  thin  cylindrical  shells  under  various 
types  of  surface  and  end  loads. 1-5 

The  investigation  presented  in  this  paper  is  concerned 
with  the  vibrations  of  a  thick-walled  visco-elastic  cylinder 
contained  by  and  bonded  to  a  thin  cylindrical  shell  (Fig.  la'-; 
the  necessity  for  considering  damping  in  the  thick-walled 
cylinder  is  due  to  the  fact  that  most  solid  propellants  ex¬ 
hibit  a  visco-elastic  response. 


The  solutions  presented  are  radially-symmetric .  Fre¬ 
quency  response  functions  are  determined  which  may  easily  be 
utilized  for  arbitrary  and  random  inputs  by  means  of  the 
established  methods  of  Harmonic  Analysis. 


Damping  is  introduced  by  means  of  the  elast ic-visco- 
elastic  analogy  and  the  use  of  complex  moduli.4  Poisson's 
ratio  v  for  the  elastic  thick-walled  cylinder  is  replaced 


by 


=  v  +  iv  and  the  shear  modulus  G  is  replaced  by 


G*  =  G  +  iG 


where 


=  G  and  G  /2ttG  denotes  the 


specific  damping  AW^/W^  ,  AW^  being  the  energy  dissipated 
during  one  cycle  in  shear  and  the  maximum  shear  strain 

energy . 

Although  the  ratio  G  G  is,  generally,  a  function  of 

2  1 

frequency,  for  the  present  considerations  it  is  assumed  to 
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be  constant  over  a  relatively  wide  range  of  frequencies. 

This  assumption  of  a  constant  complex  modulus  is  in  fair 
agreement  with  test  results.0-7 

The  dependence  of  the  damping  upon  the  amplitude  of  the 
applied  stress  is  neglected  because  its  consideration  would 
greatly  complicate  the  equations  of  the  problem  by  making 
them  nonlinear. 

The  equations  of  linear  elasticity  are  used  for  the 
thick-walled  cylinder.  Membrane  theory  is  used  for  the 
shell.  The  equations  of  a  higher  order  shell  theory  which 
includes  the  effects  of  rotatory  inertia  and  shear  are  pre¬ 
sented  in  Appendix  A.  These  equations  should  be  used  in 
cases  where  membrane  theory  is  inadequate,  e.g.,  for  external 
loads  of  acoustic  nature  or  axial  loads  distributed  over  the 
ends  of  the  shell.  Appendix  B  gives  a  solution  for  such 
axial  loads  on  the  ends  of  the  simply-supported  shell  which 
can  easily  be  superimposed  on  the  solution  for  the  shell 
under  internal  pressure. 

The  analysis  is  broken  down  into  the  following  steps: 

1.  The  thick-walled  cylinder 

2.  The  thin  shell. 

3.  The  system  of  the  thjck-walled  cylinder  contained 

by  and  bonded  to  the  thin  shell. 

The  assumption  that  the  propellant  and  case  form  a  con¬ 
tinuous  structure  at  their  interface  requires  continuity  of 
displacements  and  therefore  produces  strong  coupling  of  the 
motions  of  the  cylinder  and  shell.  This  coupling  is  con¬ 
sidered  to  be  of  primary  importance  and  therefore  the  inter¬ 
action  is  treated  rigorously. 
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As  a  first  step,  a  solution  is  obtained  for  a  thick- 
walled  visco-elastic  cylinder  with  simply-supported  ends  and 
oscillating  pressures  on  the  inner  and  outer  surfaces  as  well 
as  a  tangential  stress  on  the  outer  surface  in  the  axial  direc¬ 
tion. 

As  a  second  step,  a  solution  is  obtained  for  the  thin 
shell  with  internal  pressure  and  a  tangential  stress  in  the 
axial  direction  on  the  inner  surface.  The  outer  surface  of 
the  shell  is  tract  ion- free. 

The  third  step  is  a  combination  of  the  solutions  given 
in  steps  one  and  two  by  imposing  the  condition  that  all  dis¬ 
placements  be  continuous  at  the  interface  and  also  requiring 
the  external  loads  on  the  cylinder  to  be  compatible  with  the 
internal  loads  on  the  shell. 

Numerical  examples  are  given  which  illustrate  some  of 
the  effects  of  varying  the  geometry  of  the  structure  and  the 
specific  damping  in  the  thick-walled  cylinder. 
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II.  THE  SIMPLY-SUPPORTED  THICK-WALLED  CVLINDER 


The  equation  of  motion  of  the  elastic  thick-walled  cyl¬ 
inder  in  cylindrical  coordinates  is 

GV‘  u  •<  V7  •  u  -  ,  u  1) 

-  l-2v  -  ' - 

where 


o 


1  ± 
r  i  r 


i 

r  ^  • 


Application  of  Helmholtz's  theorem  to  the  displacement  vector 
u  allows  it  to  be  written  in  terms  of  potentials  as 

u  VI  f  v  ,  H  .  (2) 


When  (2)  is  substituted  into  1)  two  uncoupled  wave  equations 
on  and  H  are  obtained  by  separation  of  variables 


where 


v 

c  V‘'H  -  H 

s  -  — 


_  2Gj  1  -  v ) 
v  "  p( 1  -  2v ) 


(3) 


(4) 


is  the  square  of  the  velocity  of  propagation  of  dilatational 
waves  in  an  infinite  medium  and 
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is  the  square  of  the  velocity  of  propagation  of  shear  waves. 

In  the  case  of  cylindrical  symmetry  the  operator  V 
becomes 


V2 


,  I  *  4 

dr2  r  )  r  -T z 7 


and  the  components  of  displacement  are: 


u  = 


u,  = 


u  = 


ol> 

.  mao 

dr 

oz 

oH 

r 

oH 
_ z 

dz 

0  r 

OH, 

>Z 

c)r 

5  0 

for  ■ 

=  0 


<6 


The  requirement  that  u  5  0  for  the  radially-symmetric 


case  can  be  met  by  assuming  Hf  =■  H  =  0  ,  and  thus  the 
governing  equations  are  the  first  of  Eqs.  (3)  and  the  second 
with  the  component  H  y  replacing  the  vector  H  .  For  sim¬ 
plicity  H  ,  will  now  be  replaced  by  H  and  the  governing 
equations  become: 


d2<t  ,  1  d<t>  ,  d2$  1  d2$ 

drr  4  r  dz  c‘  5? 

v 

(7) 

d2H  1  dH  ,  d2H  1  d2H 

dr2  r  dr  dz2  dtr 
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The  displacements  of  interest  are  now 


_  oil 

or  oz 


,  &  4  ii 

oz  <r  r 


Use  of  these  displacements  results  in  the  stress-displacement 
relations 


’u  / u  -u  \ 

2G  r  /  r  zl 

T  - — 7-  l-\  -  \  vl —  t  - 

rr  l-2v  -r  \  r  z  / 


2G  .  \  (  ,  PLz  .  lUi 

l-2v  r  '  \  'z  '  'r  , 


’U  /  »u  u 

T  1-V)  I  v  U  4  — 

z  z  1  -  2  v  *z  l  o  r  r 


/ou  'U 

t  g(-' — —  ♦  — - 

rz  \  >z  >r 


Solutions  are  required  to  satisfy  the  following  conditions: 
On  the  ends  z  4L  2 


On  the  surface  r 


T  =0 
rz 
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and  on  the  surface  r  =  b 


Trr  ’  'Pb 


rrz  H) 


Solutions  of  Eqs .  7  are  chosen  of  the  form 


<J>(r,z,t  =  >  :■  r  cos  X^z  e 

m~ 1 


H(r,z, t  =  >  H  r  sin  \  z  e 

,  m  m 

m~l 


Substitution  of  Eqs.  (13)  into  (7)  gives 


<1 "  (  r )  4  -  $ 1  r 


-  k2  *  (r)  0 

m  m 


h"( r)  4  -  H1  (r) 

m  '  v  m  ' 


-  h  H  (r)  -  0 
m  m 


with 


The  solutions  of  either  of  Eqs.  (14)  can  be  one  of  three  types, 
depending  upon  whether  k^  or  is  negative,  zero  or 

positive.  Taking  the  equation  on  as  an  illustration, 


"7" 
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M  ( m , r )  =  2G  h  A  J  ( h  r) 

3  m  m  l  m 

H4(m,r)  =  2G  hmVl(hmr) 

Ni(ra,r)  =  2G  VmJl(kmr) 

N  <m,r)  =  2G  k  A  Y,(k  r) 

2  m  m  iv  m 

N  m,r  =  G  X‘  -  “V  -  h?)j  h  r) 

3  m  r  m  o  m 

N  (m,r)  =  G(  X'  -  -K  -  h2)Y  (h  r) 

4  m  r  mom 


(22) 


The  loads  on  the  surfaces  of  the  cylinder  are  now 
assumed  to  be  as  follows: 

p  z,t  =  P  (z)elat  =  £  P  cos  X  z  elj't 
*a  ’  1  a  am  m 

m 

Pb<z’t)  =  pb(z)ela>t  =  Pbm  cos  V  °1"t  (23) 

m 

^(z.t)  -  rh(z)eiat  -  £  Tbm  sin  V  eiu)t 

m 


where 


P 

am 


Pbm 


Tbm 


2  L/2 

L  J/2  Pa<Z)COS  V  dz 
2  L/2 

—  f  P,i  z  )cos  X  z  dz 

L  J  bv  '  m 

-L/2 
2  L/2 

—  f  T,  (z)sin  A  z  dz  . 
L  -L/2  b 


(24) 
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The  boundary  conditions  (11,  12)  now  become 


t  (a,z,t)  =  -X  Pa  cos  \  z  ela)t 
rr  am  m 

m 

t  (b.z.t)  =  -X  P,  cos  \  z  e1^ 
rr  bm  m 

m 

1 25) 


tr2(a,z,t)  =  0 

t  (b,z.t)  =  X  T  sin  >  z  ela>L  . 
rz  Dm  m 

m 

Introduction  of  the  above  conditions  25)  into  Eqs.  (21)  re¬ 
sults  in  the  system  of  simultaneous  equations  shown  in  Table 


I. 


TABLE  1 


A 

m 

B 

m 

C 

m 

D 

m 

M  (m,a) 

M  m,  a ) 

M  ‘  m,a 

M  ( m ,  a  ) 

-Pa 

l 

2 

3 

am 

M  (m,b) 

M  (  m , b) 

2 

M  (m,b) 

3 

M  m,b) 

4 

~Pbm 

N  ( m,  a ) 

N  (m,a) 

N  ( m, a ) 

3 

N  (m,a) 

4 

0 

N  (m,b) 

N  (m,b) 

2 

N  ( m,b ) 

3 

N  m,b) 

4 

T, 

bm 

The  solutions 


Tbm 


for  A  ,  B  . 
nr  m 

M  (m, a ) 

2 

M  ( m,b ) 

2 

N  ( m, a ) 

2 

N  (m_,b) 

2 


C  and  D  are  as  follows: 


m 

m 

M  (m,a) 

3 

M  (m,a) 

4 

M  (m,b) 

3 

M  (m,b) 

4 

N  (m,a) 

3 

N  (m,a) 

4 

N  (m,b) 

3 

N  (m,b) 

4 
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-P 

am 

M  (m,a 

3 

M  m,a) 

B  =  — 

M  (m,b) 

"Pbm 

M  ( m ,  b 

3 

M  (m,b‘ 

m  A 

N  (m,a) 

l 

0 

N  (m,a 

3 

N  (m,a) 

4 

N1(m,b) 

T, 

bm 

N  (m,b 

3 

N  m,b)  | 

4 

M  (m,a) 

l  ' 

M  ( m, a 

Pam 

M  m,a) 

4 

C  _  1 

M  m,b) 

M  ( m,b 

2 

_Pbm 

M  m,b) 

4 

m  A 

N  ( m,  a 

l 

N  m,a 

2 

0 

N  m,a ) 

1 

Ni(m,b) 

N  (m,b 

2 

T, 

bm 

N  (m,b) 

4 

M  (m,a 

l  ' 

M  (m,a) 

2 

M  m,a 

3 

-p 

am 

D  =  I 

Mi(m,b) 

M  m,b 

2  ' 

M  m,b ) 

3 

-P, 

bm 

m  A 

Ni(m,a) 

N  (m,a) 

2 

N  (m, a ) 

0 

Ni(m,b) 

N  (m,b) 

2 

N  (m,b) 

Tbm  . 

M  (m,a 1 

l 

M  (m, a) 

2 

M  (m,a) 

3 

M  (m,  a ) 

4 

Mi(m,b) 

M  (m,b) 

2 

M  (m,b) 

3 

M  (m,b) 

4 

N  (m,a) 

l 

N  (m,a) 

N  (m,a) 

3 

N  ( m, a ) 

4 

Ni(m,b) 

N  (m,b) 

2 

N  (m,b) 

N  (m,b) 

4 

(27) 
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It  is  also  possible  to  write  the  solutions  in  the  form 


A  P  4  A  P.  4  A  T, 
mi  am  m2  bm  ma  bm 


B  P  4  B  P.  4  B  T, 
mi  am  m2  bm  mo  bm 


CP  4  c  P  4CT 
mi  am  m2  bm  mo  bm 


D  P  4  D  P,  4  D  T, 
mi  am  m2  bm  mo  bm 


where  each  of  the  A  .  ,  B  .  ,  C  and  D  .  j  =  1,2,3)  can 

mj  mj  mj  mj  J 

be  determined  by  expansion  and  grouping  terms  in  Eqs.  ,26'. 

The  displacements  and  stresses  can  now  be  written  in 

terms  of  known  quantities  P  ,  P.  ,  T,  )  by  substitution 

n  am’  bm  bm  1 

of  Eqs.  (28)  into  (20)  and  then  using  the  results  in  Eqs.  (9) 
Thus  the  displacements  are: 

u  =  -X  <  [A  P  4-  A  P,  4  A  T.  ]k  J  ( k  r )  4 

r  1  mi  am  m2  bm  mi  bm '  m  1  m 

m  ^ 

4- [B  ,P  4  B  ,P,  4  B  T,  ]k  Y  (k  r  + 

1  mi  am  m2  bm  m »  bm  m  1  m 

+  [C  ,P  4  C  P.  4  C  T.  ]  \  J  hr  4 

L  mi  am  m2  bm  mo  bm  m  o  m 

4- [D  P  4-  D  P,  4  D  T,  ]  \  Y  h  r  l  cos  \  zeiu  l 

mi  am  m2  bm  mo  bm  m  o  m  '  (  m 

r  (29) 

u  =  -X  <  [A  P  +  A  P,  4  A  T,  |  A  J  (k  rU 

z  )  mi  am  m2  bm  mo  bm  m  o  m 

m  L 


+  [B  P  4  B  R  4  B  T,  ]A  Y  f  k  r)  4 

L  mi  am  m2  bm  mo  bm  m  o  m  ' 

4  [C  P  4  C  P,  4  C  T,  ][h  J.(h  r )  -  —  J  (h  r)]  4- 

1  mi  am  m2  bm  m3  bmJL  m  1'  m  '  r  o'  m  ' 

4  [D  P  4  D  P,  4  D  T,  ]  [h  Y.(h  r)--Y  (h  r)]lsin  Azem)t 

L  mi  am  m2  bm  m3  bmJ  L  m  iv  m  '  r  o'  m  '  (  rt 
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and  the  stresses  are 


"t  =Z  7~oT7  i  ( (  v  ( ^  —  A2 )— r )  4  “ ~\J  ( k  r  )  ]  [  A  P  4A  P.  4A  T  j  + 

rr  l-2v  1  m  m  mom  rim  mi  am  m2  bm  m3  bm 

k 

4[(v(k  -A2)-k  )Y  Ik  r)4  -^Y ,(k  r)][B  P  4B  Pv  +B  T.  ]+ 

m  m  m  o  m  r  1  mr  mi  am  m2  bm  m3  bm 

t  l-2v  h  A  J  (h  r)[C  ,P  4  C  p,  4  C  T.  ]  4 
'  m  m  1  m  mi  am  m2  bm  m3  bm 1 

4  Y.  (h  r )  (D  P  4D  P,  4 D  T,  >  cos  \  z  eiu)t 
1  m  mi  am  m2  bm  m3  bm  |  m 

T  =Z  2G  <  k  A  J,  (k  r)  [A  P  -4  A  P,  4  A  T,  1  4 

rz  1  m  m  1  m  mi  am  m2  bm  m3  bm 

m  V. 

4  k  A  Y,  (k  r )  [B  P  4B  P.  4B  T.  4 

m  m  1  m  mi  am  ms  bm  ma  bm 

t  ~  (A?-t4  h2;j  (h  r)[C  P  4  C  P,  4  c  T.  ]  4 
2  m  r  m7  o  m  7 1  mi  am  m2  bm  m3  bmJ 

4  i(  A2-i-4h2)Y  (h  r)[D  ,P  4D  Pv  4D  T,  ])sin  A  zela)t 
2  m  r  nr  o'  m  '  1  mi  am  m2  bm  m3  bmJj  m 


Tn^=-Z2G  I  [-^.(k  r)4,— A24k2)j  k  r)][A  P  4A  P,  4A  T.  ]4 
OQ  1  1  r  1  m  7  l-2v  m  m  o  m  1  mi  am  m2  bm  m3  bm 

m  ^ 

k 

4  [-^Y  ( k  r)4  T~~r  (  X24  k2  )Y  (k  r)  ]  ( B  ,P  4B  .  P.  4B  T,  ]4 
r  1  m  7  l-2v  m  m  o  m  7  mi  am  m2  bm  m3  bm 

A 

—  J  (h  r  [C  P  4  C  P,  4  C  T,  ]  4 
r  o'  m  '  mi  am  m2  bm  ma  bm 

4  —  Y  (h  r)[D  ,P  4D  P,  4D  T.  A  cos  A  z  e1U* 
r  o'  m  7  1  mi  am  m?  bm  m3  bmJ  (  m 

T  =Z2G  <  A2  -  — r~  k2 )  J  ( k  r )  [  A  P  4  A  P,  4  A  T.]4 

zz  1  Vl-2v  m  l-2v  m7  o'  m  mi  am  m2  bm  m3  bmJ 

m  V. 

+  (t^A2--~ T  k?)Y  (k  r )  [ B  ,P  4  B  P.  4  B  T.  ]4 
'l-2v  m  l-2v  m  o'  m  1  mi  am  m2  bm  m3  bmJ 

A 

4[-JBJ  (h  r)-h  A  J  (h  r)  1  [C  P  4  C  P,  4  C  T,  ]  4 

Lr  o  m  7  mmi'm/JLmiam  m2  bm  m3  bm 

A  . 

4  [— (h  r)-h  A  Y.(h  r)][D  P  4D  P.  4D  T,  ]  cos  A  ze 
1  r  o'  m  7  m  m  1'  m  '  mi  am  m2  bm  m3  bm  m 
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It  is  also  observed  that 


(m  =  1,2,3,-..) 


is  sufficient  to  satisfy  the  end  conditions  Eqs.  (10). 

Having  the  solution  for  the  elastic  thick-walled  cyl¬ 
inder  the  conversion  is  made  to  the  visco-elastic  thick-walled 
cylinder  by  replacing  the  elastic  shear  modulus  G  by  the 
complex  modulus  G*  and  Poisson's  ratio  v  by  a  complex  v* 
in  the  expressions  for  displacements  and  stresses. 

It  should  be  noticed  acre  that  if  in  Eqs.  (29)  and  (30; 

the  conditions  >  =  0  .  T.  =0  and  C  *  D  5  0  are  in- 
m  bm  mm 

troduced  the  displacements  and  stresses  become 

ur  =  [(AiPa  +  A,  pb)Ji(kr)  1  <Bipa  4  B2Pb)Y, (kr) ) le1  4 


u  2  0 
z 


Trr  =  -2G  ([AiPa  +  A.Pb1(2^I  k  Jo(kr)  +  r  Ji<kr>1  + 

+  lBipa  H  B2pbK^ri  k  Vkr)  4  r  Yi(kr)  l)  e1^ 

1ee  =  -2G  {t\pa  4  A-pb HsTrr  k  J0<kr>  ‘  r  Ji(kr)1  4 

+  lBiPa  4  B2Pb1[2^I  k  Yo(kr>  -  r  Yi<kr>)]  elmt 


-  *  f^T  {[Aipa  4  A=pb)Jo<kr)  4  'BiPa  4 
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with 


k2  =  to2/c2  and  c2  = 

It  can  then  be  shown  that  the  above  stresses 
the  equation 


G 

P  * 
(32) 


satisry 


zz 


v(T  +  T  n  i 

rr  t iO ' 


which  together  with  the  condition 


are  the  requirements  for  a  state  of  plane  strain.  Equations 
(31)  and  (32)  are  therefore  the  solution  to  the  problem  of 
the  thick-walled  cylinder  in  plane  strain. 
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III.  THE  SIMPLY-SUPPORTED  THIN  CYLINDRICAL  SHELL 


The  equations  of  membrane  theory  are  used  for  the  thin 
cylindrical  shell.3  For  solutions  which  are  radially  sym¬ 
metric  the  two  components  of  displacement  are  w  and  u  as 
shown  in  Fig.  lb. 

The  stress-displacement  relations  are 


o 


o 


_  JG_  ,w  v  du 

09  -  1-v  R  oz 


zz 


2G  /—  w  ,  Sui 

“  T^r  <v  s  +  57> 


(33) 


The  equations  of  equilibrium  are  therefore 


,  °09 

p  -  d  — — —  =  Mw 

*r  R 


da 


zz 


p  id  — \ —  =  Mii 
rz  dz 


(34) 


and  by  substituting  Eqs.  (33)  into  (34)  the  displacement 
equations  of  motion  become 


Pv  " 


2Gd  w 
(l-v)R  r 


—  du. 

v  vr)  = 

dz ' 


Mw 


2Gd  v  dw 
1-V  R  dz 


=  Mil 


(35) 


with  M  the  mass  per  unit  area  of  the  middle  surface  and  v 
and  G  Poisson’s  ratio  and  the  shear  modulus  of  the  shell. 
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It  is  now  assumed  that  the  external  loads  p  and  p 

IT  2 

are  such  that  they  can  be  represented  as  follows: 


Pr  = 

X 

m 

P,  cos 
bm 

>- 

3 

N 

eiujt 

pz  = 

-  >; 

m 

Tbm  3in 

A  z 
m 

eici)t 

(36) 


with 


,  (  2m- 1  ~  7T  .  ,  _  _ 

m  L  m  -  1*2,3*  *  * • J  • 


The  displacements  w  and  u  are  therefore  taken  to  be 
of  the  form 


u=XU  sin  A  z  e 
m  m 

m 


w  =  ^  W  cos  A  z  e 
m  m 

m 


iu.'t 


icut 


(37) 


Substitution  of  Eqs.  (36)  and  (37)  into  (35)  yields  the  sys¬ 
tem  shown  in  Table  II. 


TABLE  II 


U 

m 

Wm 

8A2  -  Mu)2 

(3  £  A 

R  m 

-T. 

bm 

~  Mcd2 

P. 

bm 
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Solutions  of  the  system  given  in  Table  II  for  u  and 

in 

are  the  following: 


U  =  A  ” 1  ( -T,  £  +  T.  Mo2  -  Pv  £  vX  ) 

m  bin  R  bm  bm  R  m 


W  =  A-1(  P,  0  X2  -  P,  mMr2  4  T,  -  "vX  ) 
m  bm  m  bm  bm  R  m 


with  (3  =  2dG/l-V  ,  d  being  the  thickness  of  the  shell,  and 
A  the  determinant  of  the  coefficients  of  the  system  given  by 
Table  II. 

If  quantities  U  ,  U  ,  W  and  W  are  now  appro- 
^  m2  m3  m2  m3 

priately  defined  from  the  coefficients  of  P^  and  T^m  in 
Eqs.  (38),  the  displacements  (Eqs.  (37))  are 


u  =  Z(U  P,  +  U  T,  ) sin  X  z  e' 
mp  bm  m3  bm  m 

m 


w  =  Z  W  P,  4  W  T,  )cos  X  z  e 
'  m2  bm  m3  bm'  m 

m 


where  P,  and  T,  are  the  Fourier  coefficients  of  the  ex- 
bm  bm 

pansions  for  p^  and  pz  as  given  by  Eqs.  (24). 

The  application  of  plane  strain  conditions 


X  2  °  ,  T,  =  0 
m  bm 


results  in  the  displacements 


u  =  0 


w  =  W  P,  e 
b 
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where 


d(2G/(l-‘v)  -  R  pa)2 

It  is  also  evident  from  the  stresses  (Eq.  33))  that 

o  =  7  (o  =0)  . 

zz  6^  '  rr  1 


(*1) 
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IV .  THE  THICK-WALLED  VISCO-ELASTIC  CYLINDER 
CASE -BONDED  TO  THE  THIN  SHELL 


Section  I)  contained  a  solution  for  the  simply-supported 
thick-walled  cylinder  with  normal  pressure  on  its  inner  and 
outer  surfaces  and  tangential  stress  on  its  outer  surface  in 
the  axial  direction.  In  Sec.  II,  a  solution  was  given  for 
the  simply-supported  shell  with  normal  pressure  and  tangen¬ 
tial  stress  in  the  axial  direction  on  its  inner  surface. 

These  two  solutions  will  now  be  combined  to  form  a  solu¬ 
tion  of  the  problem  of  the  cylinder  contained  by  the  shell. 

The  condition  that  displacements  be  continuous  at  the  inter¬ 
face  is  expressed  by  the  following  relations: 

u  (b,z,t  =  w(z,t) 

r  (42) 

uz(b,z,t)  =  u(z,t  . 

Substitution  from  Eqs.  (29)  and  (39)  into  (42)  results  in  two 

simultaneous  equations  on  Pbm  and  Tbm  which  allows  their 

determination  in  terms  of  P_  .  The  equations  are 

am 


P7  t  +  T,  t  =  -P  t, 
bm  2  bm  3  am  i 


P,  t  +  T.  t  =  -P  r. 
bm  2  bm  3  am  1 


(43) 
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with 


-22- 


The  ratio  Pfa/Pa  for  plane  strain  is  found,  either  by 
applying  the  first  Eq.  (42)  to  satisfy  continuity  of  dis¬ 
placement  using  (31)  and  (40),  or  by  imposing  the  conditions 


\ 

m 


=  0  , 


D 


m 


=  0 


on  the  first  Eq.  (46) . 
The  result  is 


Pb/Pa  =  <k  lJ0(kb)Y1(kb)  -  Yo(kb)Ji(kb)  ]  (48) 


where 

D  =  [(l-v)k  JQ(ka)  (2v-l)i  Jjfka)  ] Y , ( kb )  - 

-  [(l-v)k  Y0(ka)  +  (2v-l)i  yt(ka)  JjJkb)  -  (H9i 

_ 2R2Q\ _ 

( 2v-l  )d(  2G/(  1-v-)  -  R2po)2 


with 

A  =  ((l-v)kJo(ka)4(2v-l)jJJ(ka)][(l-v)kYo(kb)4(2v-l)gY1(kb)]  - 

-  [(l-v)kJo(kb)4(2v-l)^J1(kb) ][(l-v)kY0(ka)4(2v-l)^Yl(ka) ] 

(50) 
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V .  NUMERICAL  EXAMPLES 


In  order  to  illustrate  the  results  of  the  analysis  in 

the  previous  sections,  certain  numerical  calculations  were 

carried  out.  It  was  assumed  that  the  internal  pressure  in 

the  thick-walled  cylinder  was  of  the  form  P  eli't  with  P 

a  a 

constant.  This  latter  assumption  is,  of  course,  an  approxi¬ 
mation,  since  it  is  known  that  a  well  defined  smooth  burning 
surface  would  not  exist  even  in  a  solid  fuel  rocket  motor 
which  initially  had  a  circular  bore.lx  Therefore  the  tran¬ 
sition  region  where  tne  burning  takes  place  is  not  being  con¬ 
sidered.  The  overall  effect  on  the  remainder  of  the  pro¬ 
pellant  is  assumed  to  be  taken  care  of  by  an  equivalent  pres¬ 
sure  distribution.  The  effect  on  the  analysis  of  the  moving 
internal  boundary  is  also  neglected.  However,  a  step  by 
step  increase  of  the  ratio  a  b  is  felt  to  be  a  reasonable 
approximation  of  this  moving  boundary. 

Computations  were  made  for  the  circumferential  stress 
t  at  z  =  0  on  the  inner  surface  r  =  a )  of  the  thick- 

U'J 

walled  cylinder  and  at  the  interface  r  =  b)  for  certain  combi¬ 
nations  of  parameters.  At  r  =  b  the  quantities  t  ,  TrZ 
and  u^  were  also  computed  for  one  combination  of  parameters. 

For  the  numerical  calculations,  the  following  dimen¬ 
sionless  parameters  were  fixed: 

1)  The  ratio  of  the  length  to  the  outer  radius  of  the 
thick-walled  cylinder  (L/b  =  5-0)  • 

2)  The  ratio  of  the  shear  modulus  of  the  shell  to  that 
of  the  cylinder  (g/G  =  104 )  . 
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3)  Poisson's  ratio  for  the  thick-walled  cylinder 
(v  -  .4)  . 

4)  Poisson's  ratio  for  the  shell  ("v  =  .3)  . 

5)  The  ratio  of  the  mass  density  of  the  cylinder  to 
that  of  the  shell  ( , ,/\  ~  .5)  . 

Certain  other  parameters  were  varied  as  indicated  below. 

l)  The  ratio  of  the  inner  to  the  outer  radius  of  the 
thick-walled  cylinder  was  given  the  values 

a  “b  =  .25,  .5.  -75,  .95 

2  The  ratio  of  the  outer  radius  of  the  thick-walled 
cylinder  to  the  thickness  of  the  shell  was  given 
values 

b/d  =  500.  750,  1000 

3)  The  specific  damping  in  the  cylinder  was  varied  by 
allowing  the  ratio  of  the  imaginary  to  the  real 
part  of  the  complex  modulus  to  assume  the  values 

G  G  =  .1,  .2,  .3  ■ 

?  1 

Each  diagram  illustrates  the  frequency  dependence  of  the 
quantity  which  is  plotted.  The  independent  variable  in  each 
case  is  a  dimensionless  frequency  defined  by  the  equation 


All  calculations  were  done  on  the  IBM  7090  Date  processing 
system. 
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VI.  DISCUSSION  AND  CONCLUSIONS 


All  computations  were  made  for  m  =  1  i.e.,  for  the 

first  barrelling  mode  of  the  structure.  This  is  associated 
with  the  first  term  of  the  expansion  for  the  internal  pres¬ 
sure.  Figure  2  shows  a  comparison  of  this  m  1  term  with 
a  sum  of  the  series  and  it  appears  that  this  first  term  is  a 
good  approximation  for  the  behavior  as  a  whole  in  the  con¬ 
sidered  range. 

An  independent  computation  of  the  natural  frequencies 
of  the  thick-walled  cylinder  shows  that  in  the  range  under  con¬ 
sideration  there  are  fewer  roots  of  the  frequency  equation  than 
there  are  natural  frequencies  of  the  combined  structure  of  the 
cylinder  and  shell  Fig. 7a.  In  each  case  the  first  natural 
frequency  oi  the  shell  in  this  mode  is  beyond  the  range  of 
interest.  It  is  therefore  evident  that  additional  resonant 
frequencies  are  encountered  in  the  range  of  interest  because 
of  the  coupling  of  the  motions  of  the  cylinder  and  shell. 
Figures  3a  and  8a  indicate  this  increase  in  the  number  of 
natural  frequencies.  Figure  3b  illustrates  the  effects  of 

very  small  damping  (G  /G  .001)  on  the  circumferential 

2  1 

stress  at  the  bore  for  a/b  .5,  b  d  1000,  L/b  =  5-0  . 

It  can  be  seen  that  all  the  resonant  peaks  are  still  clearly 
defined.  Figure  3C  indicates  the  effects  of  higher  damping 
(G  /G^  =  .1,  .2,  and  .3)  on  the  same  stress  for  the  same 
geometry.  Here  it  is  clear  that  some  of  the  resonant  peaks 
have  already  been  completely  damped  out  for  the  value 

G  /G  =  .1  .  It  seems  apparent  therefore  that  since  the 

2'  1 

cylinder  is  in  its  first  breathing  mode,  the  resonant  peaks 
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which  arc  so  easily  damped  out  must  be  associated  with  sec¬ 
ondary  shearing  effects  that  are  coupled  with  the  barrelling 
motion.  For  higher  values  of  damping  G  'G  =  .2  and  .  j 
there  still  remains  at  leastone  prominent  resonant  peak. 

This  is  associated  with  the  first  resonant  frequency  of  the 
breathing  mode  of  the  structure. 

Figures  7  to  l6a,b,c  illustrate  the  effect  of  changing 
geometry  on  the  response  in  the  range  of  interest.  It  is 
clear  that  as  a  "b  •  1  Figs.  3C <  7b.  bb.  ;)b  the  dominant 
resonant  frequency  shifts  towards  a  value  which  is  the  first 
resonant  frequency  of  the  shell  by  itself.  As  the  change 
takes  place,  i.e.,  as  the  visco-elastic  cylinder  gets  rela¬ 
tively  thinner  it  is  also  clear  that  the  amplitude  of  the 
response  at  the  resonant  frequency  increases  in  magnitude 
tending  towards  the  limiting  case  of  the  infinite  value  which 
would  be  associated  with  the  elastic  shell  at  its  natural 
frequency.  A  similar  shifting  of  the  frequencies  can  be  ob¬ 
served  as  b  d  varies  from  1000  to  500  Figs.  7b,  10b,  l4b). 
This  amounts  to  a  relative  thickening  of  the  shell.  In  this 
case  however,  the  effects  of  damping  do  not  disappear  as 
fast,  since  in  this  case,  the  relative  thickness  of  the 
cylinder  is  held  constant  while  the  shell  is  made  thicker. 

A  study  of  Fig.  7  to  l6b  would  therefore  indicate  that  there 
is  a  certain  set  of  ratios  a  T)  and  b  d  for  which  the 
maximum  effect  of  the  damping  in  the  thick-walled  cylinder 
would  be  achieved. 

Figures  jd  and  7  to  lbc  illustrate  the  effects  of  damp¬ 
ing  on  the  circumferential  stress  at  the  interface  (r  =  b). 

It  is  clearly  noticeable  that  in  certain  cases  (Figs.  3d, 

10c  and  l4c )  where  the  cylinder  is  still  much  thicker  than 
the  shell,  the  amplitude  of  the  stress  at  the  interface  is 
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very  much  smaller  than  at  the  bore  even  though  the  resonant  peaks  are 
still  cvi cent.  however, as  is  seen  in  Figs,  jc ,  12c  and  l6c  as 
a/h  -*  1  the  amplitude  again  becomes  very  large  as  the  reso¬ 
nant  frequency  approaches  that  of  the  elastic  shell  by  itself. 

Figures  4a, b  and  oa.b  show  one  case  of  the  radial  stress 

and  displacement  at  the  interface  for  z  0.  while  Fig.  |>a, b 

shows  the  shear  stress  at  the  interface  for  z  i L  2  .  As 

would  be  expected  for  this  mode  the  shear  stress  is  much 

smaller  in  amplitude  than  the  circumferential  or  radial 

stresses.  However,  it  is  interesting  to  note  that  in  this 

case  there  is  still  evidence  of  all  the  peaks  for  G  G  .1 

2  1 

whereas  for  the  other  stresses  most  of  these  have  already 
been  completely  damped  out.  Therefore,  there  seems  to  be 
evidence  here  to  support  the  earlier  conclusion  that  those 
peaks  in  the  other  stresses  which  disappeared  for  very  low 
damping  were  due  to  frequencies  associated  with  shearing. 

It  is  apparent  from  the  results  discussed  above  that 
values  of  G  G  ^  .1,  .2,  .3  have  a  drastic  effect  upon 

the  amplitudes  of  tne  frequency  response  functions  which 
were  computed.  Since  these  values  are  smaller  than  values 
measured  experimentally  G  G  .oiu)1  it  seems  reasonable 
to  conclude  that  any  analysis  of  a  solid  fuel  rocket  motor 
which  does  not  take  into  consideration  the  visco-elastic 
effects  of  the  propellant  would  not  give  a  true  overall  pic¬ 
ture  . 

The  variation  of  the  ratio  a  "b  ,  as  a  simulation  of 
the  state  of  the  structure  at  various  stages  of  burning, 
indicates  for  certain  ranges  of  frequency,  a  considerable 
increase  in  the  amplitudes  of  the  computed  functions.  This 
implies  that,  in  fact,  in  a  structure  of  this  type,  the 
critical  period  would  be  in  the  later  stages  when  the  pro- 


pellant  is  almost  completely  burnt  out.  From  this  point  of 
view,  the  important  frequency  might  be  close  to  the  first 
breathing  frequency  of  the  thin  shell. 

However,  the  diagrams  indicate  that  for  certain  re¬ 
stricted  ranges  of  the  forcing  frequency  a  judicious  choice 
of  parameters  may  be  enough  to  eliminate  the  possibility  of 
large  amplitude  resonant  peaks. 
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APPENDIX  A 


H IGHER  ORDER  SHELL  THEORY 


In  cases  of  dimensions  and  loads  for  which  membrane 
theory  is  inadequate,  a  higher  order  theory  is  required.  The 
theory  presented  here  includes  the  effects  of  rotatory 
inertia  and  shear.13 

For  the  radially  symmetric  case  the  displacements  are 


U  (x,2,t)  =  W(z,t) 

r  (Al) 

uz(x,z,t  =  u(z,t  +  X^(z,t' 


with  u,  w,  as  shown  in  Fig.  lb.  x  =  r-R  and  V'  the  angle 
of  rotation  of  fibers  which  in  the  initial  state  are  perpen¬ 
dicular  to  the  middle  surface. 

The  relevant  strains  are 


f 


00 


€ 

zz 


rz 


,  du 
-  2- 
“  2  Sx 


au 

r 

az 


(A2) 


Although  unrestricted  Poisson's  ratio  effect  is  allowed 
in  the  radial  direction  and  at  the  same  time  the  displace¬ 
ments  seem  to  imply  =  0  >  the  contradiction  can  be  re¬ 

moved  by  first  eliminating  e  from  the  stress-strain  re¬ 
lations  . 
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Thus 


with 


(A5) 


p 


2dG 

"  i-v  ; 


d^G 
6( i-v' 


K?dG  . 


All  other  forces  and  moments  are  zero.  Application  of  the 
necessary  conditions  for  equilibrium  of  a  section  of  the 
shell  results  in  three  uncoupled  equations  on  u,  w,  and  y 


v  2  ^ 

P(R  +  v  ^)  +  R  pz  -  RMli  =  0 

4  v  fl>  -  »(R  4  R  0)  -  R  Px  '  ■ 0  <A6> 

A  a  '  /  3  w  d  •  • 

6R  £  -  7(Ry  +  R  §7)  ■  R  2  PZ  ■  RIv  "  0 

where  M  is  the  mass  and  I  the  moment  of  Inertia  per  unit 
area  of  the  middle  surface  of  the  shell. 

It  can  be  shown  that  a  unique  solution  of  the  system 
(A6)  can  be  obtained  if  one  quantity  of  each  of  the  pairs 
Nzu,  Mz V7  and  Qzw  is  prescribed  on  boundaries  z  =  constant. 

It  can  also  be  shown  that  Eqs.  (35  can  obtained  very 
easily  from  (a6)  by  letting  y  =  -^w/t3z  and  neglecting  ro¬ 
tatory  inertia  (I  =  0)  and  bending  rigidity  (6  =0).  If  the 
loading  is  the  same  as  given  by  (36)  in  Sec.  II,  the  solu¬ 
tion  procedure  is  analogous,  the  only  difference  being  that 
in  the  case  of  the  higher  order  theory,  the  expressions  for 
tj  in  (44)  will  contain  an  additional  term  due  to  ^  . 
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Numerical  calculations  were  made  using  the  higher  order 
theory  in  this  section  for  the  same  range  as  for  the  simple 
shell  theory,  and  it  was  found  that  in  this  range,  the  re¬ 
sults  were  identical.  This  agrees  with  the  results  of 
Greenspon. lc 
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APPENDIX  B 

SOLUTIONS  OF  HIGHER  ORDER  SHELL  EQUATIONS 
FOR  LONGITUDINAL  LOADS  DISTRIBUTED  OVER  THE  ENDS 


Solutions  to  Eqs .  A6)  are  taken  of  the  form 


u  =  Ue 


w  =  We 


i  (  XZ+CUt  ) 
1  (  \Z+(JL>t 


i  e 


i  (  Xz+cut ) 


Substitution  of  Eqs.  (Bl)  into  Eqs.  (A6)  gives  a  system  of 
equations  on  U,  W  and  'F  which  is  homogeneous  Px  =  0  , 


=  0)  as  shown  in  Table  B-i  , 


TABLE  B-I 


In  order  to  have  a  non-trivial  solution  of  the  system 
given  in  Table  B-I,  the  requirement  is  that  the  determinant 
of  the  coefficients  of  U,  W  and  i  be  equal  to  zero. 

This  expanded  determinant  is  the  frequency-wavelength  rela¬ 
tionship  for  the  shell. 

If  a  dimensionless  wavenumber  \  =  d'  is  defined  it  can 
be  shown  that  there  are  three  limiting  values  of  the  fre¬ 
quency  u5  as  this  wavelength  approaches  zero.  These  values 
are  determined  by  the  roots  of  the  equation: 


S5° 


‘  ^  ■*  1-v")  l“  + 


2tt  *  1-v 


or  -  0 


(B2 


with 


d"  -  dVR* 


and 


tu?d  ’P 

<n  “  iJFcr 


The  roots  in  ascending  order  of  magnitude  are 


(JD‘ 

i 


2n  '1-v ■ 


71?  =  4r 

,  ir 


(B3 


The  first  root  -  0/  corresponds  to  a  rigid  body 

displacement  of  the  shell  in  the  axial  direction.  The  re¬ 
mainder  of  the  branch  which  passes  through  this  point  gives 
frequency-wavelength  relations  for  purely  axial  motions  ex¬ 
tensions)  of  the  middle  surface  of  the  shell. 

The  second  root  (u/5  =  dr'/27T2 ( 1-v )  )  corresponds  to  the 

e 

lowest  "breathing"  mode  of  the  shell.  The  motion  in  this 
mode  can  be  compared  to  the  vibrations  of  a  circular  ring  in 
its  plane.18 
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The  equation  for  the  actual  frequency  is  in  this  case 


co 


2 

2 


2 '  1+ v  ,  G 
=  R*‘l-Vz) 


(  b4  ) 


where  refers  to  Timoshenko's  value  for  the  frequency  of 

the  lowest  mode  of  a  vibrating  ring.lb  The  factor  (1-v2) 
appearing  in  the  denominator  is  a  correction  which  is  due  to 
the  effect  of  the  adjacent  rings  of  the  shell  upon  each  other. 

It  should  also  be  noted  that  -tim  co2  «=  0  indicating 

00  * 

plate . 

The  third  root  (a?  =  3/7r?)  corresponds  to  the  cut-off 

3 

frequency  for  axial  shear  in  the  shell. 

The  expansion  of  the  determinant  given  in  Table  B-I, 
has  three  roots  V  which  are  frequency  dependent. 

In  the  range  where  these  are  real,  the  solutions  (B-I 
become 

u  =  l  Uv  ei(±\(a))z+u* 
k=l  k 

w  =  X  w,  e1  -  k  '  z+  ■t  ; 35) 

k=l 


t  =  I  ei(  +  Ak(aOz+a)t) 

k-l  k 


It  is  also  possible  to  determine  ratios 


uk  3viAkM 

=  RpXz(cjo)  -  RMo: 
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(b6) 


\  to) 

Wk  IuT  -  7  -  6\£(ui) 

Symmetric  portions  (symmetry  with  respect  to  z  =  0  of  the 
solution  are  therefore 


u 


3  ,»v  sin  \  a))2 

k:L  T  -  t^k(a>; 


iu)t 

e 


w 


X  W.  cos  \  .u )  eijt 
k-1  K  K 


v  - 


y\  “)wk sln  Vu  z  iu>t 


A  - - TTTT - C 

_  ^  7  -  1(0  -  0  ( O) ) 


(B7) 


Using  these  latter  equations  it  is  now  possible  to  determine 
eigenfunctions  and  natural  frequencies  for  symmetric  motions, 
by  imposing  the  necessary  boundary  conditions. 

Using  Eqs.  (B7,  above  and  Eqs.  A5  the  following 
boundary  conditions  are  applied: 


w  =  0 


Mz  =  0  \ z  =  4L/2 


N 


z 


N  e 


icotj 


(b8) 


Substitution  of  Eqs.  ( B7 )  into  (b8)  results  in  the  system 
shown  in  Table  B-II. 
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TABLE  B-II 


This  system  can  now  be  solved  for  k  =  1,2,3)  in 

terms  of  M  the  amplitude  of  the  applied  load  as  follows: 


M  f  7*“  I 

W  =  -  { — — 7  g  ,  ,-C2'  -  T~  i  ,  c  .1>  /  COS  A  L  2  cos  A  L  '2 

i  A  I  7  -  Iur  +  6  A  7  -  Iur  +  0  A  2  u 

^  2  3> 


N  \  *VA  7*T  1 

w  =7-  { - T  2  ,  ■■  ^7 - T  '  ^  V  11?  /cos  a  L/2  cos  X  L/2  By) 

2  A  ^7  -  Iur  +  6A  7  _  Icu  +  '  3  1 


r  7^2  7a^  'i 

W  =7  - T  j  rvr - T  h  ,  V\P'  >  cos  A  L  2  cos  A  L  '2  . 

3  A  1  7  -  lor  +  6  A  7  -  Iur  +  0  A"  [  1  2 


These  quantities  can  then  be  substituted  into  Eqs .  (B7)  and 
then  combined  by  Eqs.  (Al)  to  give  the  displacements  in  the 
shell  as 
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u  =  I  W,  cos 
r  k-i  k 


Xk(!i>  z  e 


irnt 


U2  =  * 
k=l 


hv 


R.Mo-  -  ,<  \*  a 


(BIO) 


x 


7  -  IlD^  4  6\^(d)) 


Xk  a  Wk  sin  k  a‘  z  /  e 


ia)t 


Since  the  problem  is  linear,  the  results  above  can  be 
superimposed  upon  those  given  m  Sec.  II  to  be  used  as  a 
solution  for  the  combined  loading  of  internal  pressure  on 
the  thick-walled  cylinder  and  longitudinal  loads  on  the  ends 
of  the  shell. 

The  super-position  could  be  carried  out  by  first  ex¬ 
panding  the  results  BIO  into  Fourier  Series  and  then  com¬ 
bining  the  results  with  the  series  given  by  Eqs.  (39)- 
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Figure  3b  Absolute  Value  of  Circumferential  Stress  at  the  Bore 
for  Small  Damping  a  b  =  .5  ,  b  d  *  1000! 


Figure  3C  Absolute  Value  of  Circumferential  Stress  at  the  Bore 
for  Larger  Damping  (a  b  =  .5  ,  b  d  *  lOOO) 


Figure  3d  Absolute  Value  of  Circumferential  Stress  at  the 
Interface  for  Various  Values  of  Damping 
a  b  =  . 5  ,  b  d  =  1000' 
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Figure  4a  Radial  Stress  at  the  Interface  for  Elastic  Cylinder 
and  Shell  (a  b  =  .5  ,  b  d  -  1000' 
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* igure  5a  Shear  Stress  at  the  Interface  for  the  Clastic  Cyl¬ 
inder  and  Shell  (ab=  .5  .  b  d  ™  1000) 
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Figure  6b  Absolute  Value  of  Radial  Displacement  at  the  Inter 
face  for  Various  Values  of  Damping  (ab  ■  .5  . 


Figure  7a  Circumferential  Stress  at  the  Bore  far  the  Elast 
Cylinder  and  Shell  and  also  Showing  Roots  of  t 
Frequency  Equation  for  the  Cylinder  Alone 


Figure  'fh  Absolute  Value  of  Circumferential  Stress  at  the  Bo: 

for  Various  Values  of  Damping  and  Indicated  Valut 


Figure  6b  Absolute  Value  of  Circumferential  Stress  at  the 
for  Various  Values  of  Damping  and  Indicated  V 


Figure  8c  Absolute  Value  of  Circumferential  Stress  at  the  In 
terface  for  Various  Values  of  Damping  and  Indi- 
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Figure  2b  Absolute  Value  of  Circumferential  Stress  at  the 

Bore  for  Various  Values  of  Damping  and  Parameters 
as  Indicated 


Figure  9c  Absolute  Value  of  Circumferential  Stress  at  the 


Figure  10a  A  Circumferential  Stress  at  the  Bore  for  the 

Elastic  Cylinder  and  Shell  and  Parameters  as 
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Figure  10b  Absolute  Value  of  Circumferential  Stress  at  the 

Bore  for  Various  Values  of  Damping  and  Parameters 


Figure  10c  Absolute  Value  of  Circumferential  Stress  at  the 
Interface  for  Various  Values  of  Damping  and 
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Figure  11a  Circumferential  Stress  at  the  Bore  for  the  Elastic 
Cylinder  and  Shell  and  Parameters  as  Indicated 


Figure  lib  Absolute  Value  of  Circumferential  Stress  at  the 

Bore  for  Various  Values  of  Damping  and  Parameters 
as  Indicated 
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Figure  12c  Absolute  Value  of  Circumferential  Stress  at  the 
Interface  for  Various  Values  of  Damping  and 
Parameters  as  Indicated 
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Figure  14a  Circumferential  Stress  at  the  Bcre  for  the  Elastic 
Cylinder  and  Shell  and  Parameters  as  Indicated 


Figure  14b  Absolute  Value  of  Circumferential  Stress  at  the 

Bore  for  Various  Values  of  Damping  and  Parameters 


Figure  14c  Absolute  Value  of  Circumferential  Stress  at  the 
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Attn:  Capt.  W.  T.  Sawyer,  Code  406 


Chiet.  Bureau  of  Ships 
Department  of  the  Navy 
Washington  25.  D.C. 

Attn:  Dr.  George  Sponsler,  Code  jl5 

Di rector 

Naval  Reserach  Laboratory 
Department  of  the  Navy 
Washington  25.  D.C. 

Attn:  Mr.  P.  Waterman.  Code  536c 

Missile  and  Space  Division 
Lockheed  Aircraft  Corporation 
Palo  Alto.  California 
Attn:  Dr.  W.  F.  Whitmore 

Special  Projects  Office  SP-11^ 
Bureau  of  Naval  Weapons 
Department  of  the  Navy 
Washington  25.  D.C. 

Attn:  LCDR  R.  H.  Yerbury 

Executive  Secretary 


PART  II  -  Members,  SPIA  Physical 
Properties  Panel 


Commander 

Air  Force  Flight  Test  Center 
Edwards  Air  Force  Base,  California 
Attn:  FTRS ,  D.  Hart 

Aberdeen  Proving  Group 
Ballistic  Research  Laboratories 
Aberdeen,  Maryland 
Attn:  A.  S.  Elder 

H.  P.  Gay 

Redstone  Arsenal 

Army  Rocket  and  Guided  Missile  Agency 
Huntsville,  Alabama 
Attn:  T.  H.  Duerr 

Department  of  the  Navy 
Bureau  of  Naval  Weapons 
Washington  25,  D.C. 

Attn:  RMMP-22,  W.  A.  Bernett 


U.S.  Naval  Ordnance  Test  Station 
China  Lake,  California 
Attn:  K.  H.  Bischel 
A.  Adlcoff 

U.S.  Naval  Propellant  Plant 
Indian  Head,  Maryland 
Attn:  W.  J.  Marciniak 

Aerojet-General  Corporation 
P.O.  Box  296 
Azusa,  California 
Attn:  K.  H.  Sweeny 
K.  W.  Bills 

Aerojet-General  corporation 
P.O.  Box  1168 
Sacramento ,  California 
Attn:  J.  H.  Wiegand 

Atlantic  Research  Corporation 
Shirley  Highway  &  Edsall  Road 
Alexandria,  Virginia 
Attn:  M.  G.  DcFries 

California  Institute  of  Technology 
Pasadena,  California 
Attn:  J.  P.  Blatz 

M.  L.  Williams 

E.  I.  duPont  de  Nemours  &  Co. 

Gibbs town,  New  Jersey 
Attn:  R.  D.  Spangler 

Grand  Central  Rocket  company 
P.O.  Box  111 
Redlands,  California 
Attn:  E.  Fitzgerald 

Hercules  Powder  Company 
Allegany  Ballistic  Laboratory 
Cumberland,  Maryland 
Attn:  J.  H.  Thacher 

Hercules  Powder  Company 

Bacchus  Works 

Magna,  Utah 

Attn:  D.  E.  Nicholson 


Jet  Propulsion  Laboratory 
4800  Oak  Grove  Drive 
Pasadena  3»  California 
Attn:  R.  F.  Landel 

Rocketdyne 

Solid  Propulsion  Operations 
P.O.  Box  548 
McGregor ,  Texas 
Attn:  S.  C.  Britton 

Rohm  and  Haas  Company 
Redstone  Arsenal  Research  Division 
Huntsville,  Alabama 
Attn:  A.  J.  Ignatowski 

Space  Technology  Laboratory.  Inc. 
5730  Arbor-Vitae  Street 
Los  Angeles  45,  California 
Attn:  w.  G.  Gottenberg 

Stanford  Research  Institute 
Menlo  Park ,  California 
Attn:  Dr.  T.  L.  Smith 

Thiokol  Chemicals  corporation 
Redstone  Division 
Huntsville,  Alabama 
Attn:  M.  H.  Cooper 

United  Technology  Corporation 
P.O.  Box  358 
Sunnyvale,  California 
Attn:  Dr.  Iwanciow 

Dr.  F.  Lavacot 

Solid  Propellant  Information  Agency 
A PL,  JHU,  8621  Georgia  Avenue 
Silver  Spring,  Maryland 
Attn:  M.  T.  Lyons 


PART  III  -  Activities  and  Contractors  -  Concern* d  with  Propellant  Mechanic;. 
Government 


Chief  of  Naval  Research 
Department  of  the  Navy 
Washington  25,  D*  C. 

At  t  , :  Code  439  3 

Code  4 11 

Commanding  Officer 
Cognizant  CNR  Bra* c)  Office 

Armed  Services  Technical 
Information  Agency 
Arlington  Hall  Station 
Arlington  12,  Virginia  (10 

Office  of  Technical  Services 
Department  of  Commerce 
Washington  25.  D.C. 

Director  of  Defense  Research  and 
Engineer j  ng 
The  Pentagon 
Washington  25,  D.C. 

Attn:  Technical  Library 

Advanced  Research  Projects  Agency 
Defense  Research  and  Engineering 
The  Pentagon 
Washington  25.  D.C. 

Attn:  A.M.  Rufcenstein 

Director .  Special  Projects 
Department  of  the  Navy 
Washington  25.  D.C. 

Attn:  SP001  (Dr.  J.  P.  Craven 

SP2J1  ( LCDR  R.  L.  McArfby 

Chief,  Bureau  of  Naval  Weapons 
Department  of  the  Navy 
Washington  25,  D.C. 

Attn:  RRRE-6  (Dr.  C.  Boyars) 

RMMP-2  vDr.  0.  H.  Johnson 
RMMP-11  Mr.  I.  Silver, 

Commander 

Air  Force  Flight  Test  Center 
Edwards  Air  Force  Base,  California 
A^tn :  FTRS 

Commander 

Ogden  Air  Material  Area 
Hill  Air  Force  Base,  Utah 
Attn:  00MQCC ,  H.  A.  Mat is 


Commander 

Air  Fore*  Office  of  Scientific  Research 
Washing  tor.  D.C. 

Art:.:  Mechanics  Division 

Commanding  General 
Aberdeen  Proving  Grou,.d 
Mary  la-  d 

Attn:  Balli:  tic  Research  Labe. 

ORDBG-BLI 

Comma -der 

Army  Rocket  and  Guided  Missile  Agency 
Red.  tor .o  Arrcnal  ,  Alabama 
Attn:  Technical  Library 
ORDXR-OTL 
ORDAB-HSI 

Department  of  the  Army 
Office.  Chief  of  Ordnance 
Washington  25.  D.C. 

Attn:  ORDTB .  J.  A.  Chalmers 

Director 

Pla: tics  Tech.  Eval.  Center 
Pi  it.  i:.ny  Ar..enul 
D  i  New  Jersey 

U.S.  Army  Research  Office 

2127  Myrtle  Drive 

Dak*  Station 

Durham  North  Carolina 

Attn:  Div.  of  Engineering  Sciences 

Chief  of  Nival  Operations 
Department  of  the  Navy 
Washington  25.  D.C. 

Atti  :  Op  07T 
Op  G3EG 

Quality  Evaluation  Laboratory 
Naval  Ammunition  Depot 
Concord,  California 
Attn:  D.  R.  Smathers 

U.S.  Nav.il  Ordnance  Laboratory 
Non  Metallic  Materials  Division 
Silver  Spring,  Maryland 
Attn:  H.  A.  Perry,  Jr. 

U.S.  Naval  Ordndnce  Test  Station 
China  Lake,  California 
Attn:  J.  T.  Bartling 


U.S.  Naval  Propellant  Plant 
Indian  Head,  Maryland 
Attn:  J.  Browning 

L.  Papier  VL) 

Office  of  Naval  Research 
Branch  Office 
495  Summer  Street 
roston  10,  Massachusetts 
Attn:  Dr.  J.  H.  Faull,  Jr. 

National  Aeronautics  &  Space  Adm. 
1570  H  Street,  N.W. 

Washington  25,  D.C. 

Attn:  Chief,  Div.  of  Research 
Information 


Contractors 


Atlantic  Research,  Inc. 

Shirley  Highway  &  Edsall  Road 
Alexandria,  Virginia 
Attn:  M.  G.  DeFries 

Battelle  Memorial  Institute 
505  King  Avenue 
Columbus  1,  Ohio 
Attn:  J.  Harry  Jackson 

Brown  University 
Division  of  Applied  Mathematics 
Providence  12,  Rhode  Island 
Attn:  Prof.  E.  H.  Lee 

Prof.  R.  S.  Rivlin 

University  of  California 
College  of  Engineering 
Berkely  4,  California 
Attn:  Prof.  Paul  M.  Naghdi 

Catholic  University  of  America 
Department  of  Civil  Engineering 
620  Michigan  Avenue,  N.E. 
Washington,  D.C. 

Attn:  Prof.  A.  J.  Durelli 

Columbia  University 
Department  of  Civil  Engineering 
Amsterdam  Avenue  &  120th  Street 
New  York  27,  New  York 
Attn:  Prof.  A.  M.  Freudenthal 


Materials  Technology,  Inc. 

11  Leon  Street 

Boston  15.  Massachusetts 

Attn:  Dr.  R.  G.  Cheatham 

New  York  University 

Depart,  of  Aeronautical  Engineering 

University  Heights 

New  York  53.  New  York 

Attn:  Prof.  H.  Becker 

University  of  Pennsylvania 
Graduate  Division  of  Engineering  Mechanics 
Philadelphia  4,  Pennsylvania 
Attn:  Prof.  Z.  Hashin 

Polytechnic  Institute  of  Brooklyn 
333  Jay  street 
Brooklyn  1,  New  York 
Attn:  Prof.  F.  Romano 

Prof.  J.  Klosner 
Prof.  F.  Ullman 

Southwest  Reseatch  Institute 
85OO  Culebra  Road 
San  Antonio  6,  Texas 
Attn:  Dr.  R.  C.  DeHart 

Central  Laboratory  T.  N.  O. 

134  Julianalaan 
Delft,  Holland 
Attn:  Dr.  F.  Schwarzl 

Aerojet-General  Corporation 
Solid  Rocket  Plant 
Sacramento.  California 

Attn:  Dr.  W.  0.  Wetmore  (2) 

Aerojet-General  Corporation 
P.0.  Box  1168 
Sacramento,  California 
Attn:  A.  Fraser 

Dr.  Zickel 

Amcel  Propulsion,  Inc. 

Box  3049 

Asheville,  North  Carolina 
Attn:  R.  N.  Lowrey 

American  Cyanamid  Company 
1937  West  Main  Street 
Stamford,  Connecticut 
Attn:  Dr.  V.  Wystrach 


Chief,  Bureau  of  Ships 
Department  of  the  Navy 
Washington  c5,  D ,  C . 

Attm  Code  335 

professor  R.  L.  Bisplinghoff 
Department  of  Aeronautical  Engineering 
Massachusetts  Institute  of  Technology 
Cambridge  39>  Massachusetts 

Chief,  Bureau  of  Yards  tc  Docks 
Department  of  the  Navy 
Washington  25,  D.  C. 

Attm  Code  70 

Commanding  Officer  x  Director 
David  Taylor  Model  Basin 
Washington  7>  C. 

Attm  Code  700 

Director 

Materials  Laboratory 
New  York  Naval  Shipyard 
Brooklyn  1,  New  York 

Of ficer- in-Charge 

Naval  Civil  Engineering  Research 

ana  Evaluation  Laboratory 

U .S .  Naval  Construction  Battalion 

Center 

Port  Hueneme,  California 
Commander 

U.S.  Naval  Proving  Ground 
Dahlgren,  Virginia 

Commanding  Officer  v  Director 
b  .S .  Naval  Engineering  Experiment 
Station 

Annapolis,  Maryland 
Super in  tendent 

U.S.  Naval  Postgraduate  School 
Monterey,  California 

National  Sciences  Foundation 
1520  H  Street,  N.W . 

Washington,  D .  C . 

Attm  Engineering  Sci.  Div . 


Professor  H.  H.  Bleich 
Department  of  Civil  Engineering 
Columbia  University 
Amsterdam  Avenue  k  120th  Street 
New  York  27,  New  York 

Professor  B.  A.  Boley 
Department  of  Civil  Engineering 
Columbia  University 
Amsterdam  Avenue  &  120th  Street 
New  York  27 ,  New  York 

Professor  B.  Budiansky 
Department  of  Mechanical  Engineering 
School  of  Applied  Sciences 
Harvard  University 
Cambridge  38,  Massachusetts 

Professor  G.  F.  Carrier 
Pierce  Hall 
Harvard  University 
Cambridge  38,  Massachusetts 

Professor  D.  C.  Drucker 
Division  of  Engineering 
Brown  University 
Providence  12,  Rhode  Island 

Professor  J.  Encksen 
Mechanical  Engineering  Department 
Johns  Hopkins  University 
Baltimore  18 ,  Maryland 

Professor  A.  C.  Eringen 
Dopciriment  of  aeronautical  Engineen 
Purdue  University 
Laf.iyette,  Indiana 

Mr.  Martin  Goland,  President 
Southwest  Research  Institute 
85JU  Culebra  Road 
San  Antonio  0,  Texas 

Professor  J.  N.  Goodier 
Department  of  Mechanical  Engineering 
Stanford  University 
Stanford,  California 


University  of  California 
Herkely,  California 
Attn:  Dr.  K.  S.  Pister 

Catholic  University  of  America 
Department  of  Civil  Engineering 
Washington,  D.c. 

Attn:  J.  Baltrukoms 

University  of  Florida 
College  of  Engineering 
Gamscvi  lie,  Florida 
Attn:  J.  Griffith 

University  of  Illinois 
Department  of  Aero  Engineering 
Urbana,  Illinois 
Attn:  Dr.  H.  H.  Hilton 

E.  I.  duPont  do  Nemours  and  Co. 
Gibbstown.  New  Jersey 
Attn:  R.  D.  Spangler 

Grand  Central  Rocket  Company 
P.0.  Box  111 
Redlands,  California 
Attn:  A.  T.  Camp 

Hercules  Powder  Company 
Allegany  Ballistics  Laboratory 
Cumberland.  Nary  lane. 

Attn:  Dr.  R.  Steieberger  ? 

Jet  Propulsion  Laboratory 
4300  Cax  Grove  Drive 
Pasadena  3.  California 
Attn:  G.  Lewis 

Lockheed  Missile  V  Space  Company 

1122  Jagcls  Road 

Sum  lie,  California 

Attn:  E.  Luken  (2) 

North  American  Aviation 
Rocketdyne  Division 
66 33  Canoga  Avenue 
Canoga  Park.  California 
Attn:  F.  Cramer 

Rohm  (c  Haas  Company 
Redstone  Arsenal  Research  Division 
Huntsville,  Alabama 
Attn:  H.  Shuey 


Thiokol  Chemical  Corporation 
Redstone  division 
Huntsville,  Alabama 
Attn:  Technical  Director 
J .  Wise 


TART  IV  -  Activities  and  contractor 
concerned  with  other  aspects  of 
Elastomer  Mechanics 


Government 

Commanding  Officer 
Office  of  Naval  Research 
Branch  Office 

John  Crernr  Library  Bui?d*ng 
o  E.  Randolph  Street 
Chicago  11.  Illinois 

Commanding  Of f icer 

Of  lice  of  Naval  Research 

branch  Office 

i» 6  Broadway 

New  York  lj.  New  York 

Commanding  Officer 
Office  of  Naval  Research 
Branch  Office 
10jJ  E.  Green  Street 
Pasadena,  California 

Commanding  Officer 
Office  of  Naval  Research 
Branch  Office 
1000  Geary  Street 
San  Francisco.  California 

Commanding  Officer 
Office  of  Naval  Research 
Branch  Office 

Navy  1  0.  Fleet  Post  Office 
Box  j j  F PO 

New  York.  New  York  (5y 

Director 

Naval  Research  Laboratory 
Washington  26.  D.C. 

Attn:  Tech.  Info.  Officer  (6' 

Code  6200 
Code  6210 


Professor  P.  G.  Hodge 
Department  of  Mechanics 
Illinois  Institute  of  Technology 
Chicago  16,  Illinois 

Professor  N.  J.  Hoff,  Head 
Division  of  Aeronautical  Engineering 
Stan  for'*  mvorsity 
Stanford,  California 

Professor  J.  Kempner 

Dept,  of  Aeronautical  Engineering 

and  Applied  Mechanics 

Polytechnic  Institute  of  Brooklyn 

333  day  Street 

Brooklyn,  New  York 

Professor  R.  D.  Mindlin 
Dept,  of  Civil  Engineering 
Columbia  University 
Amsterdam  Avenue  fc  120  Street 
New  York  27,  New  York 

Professor  William  A.  Nash 
Dept,  of  Engineering  Mechanics 
University  of  Florida 
Gainesville,  Florida 

Professor  N.  M.  Ncwmark,  Head 
Dept,  of  Civil  Engineering 
University  of  Illinois 
Urbana,  Illinois 

Professor  E.  Reiss 

Institute  of  Mathematical  Sciences 
New  York  University 
25  Waverly  Place 
New  York  3,  New  York 

Professor  W,  Prager,  Chairman 
Physical  Sciences  Council 
Brown  University 
Providence  12,  Rhode  Island 

Professor  E.  Reissner 
Dept,  of  Mathematics 

Massachusetts  Institute  of  Technology 
Cambridge  39,  Massachusetts 

Professor  Bernard  W.  Shaffer 
Dept,  of  Mechanical  Engineering 
New  York  University 
University  Heights 
New  York  53,  New  York 


Professor  Eli  Sternberg 
Dept,  of  Mechanics 
Brown  University 
Providence  12,  Rhode  Island 

Professor  A.  S.  Velestos 
Dept,  of  Civil  Engineering 
University  of  Illinois 
Urbana,  Illinois 


